
alpha(1:3)=sqrt((103+256+31)/(103+31+256*1));

alpha(4)=sqrt((103+256+31)/(103+31+256*0.783493287693));

alpha(5)=sqrt((103+256+31)/(103+31+256*0.650364580476000));

alpha(6)=sqrt((103+256+31)/(103+31+256*0.581489562953000));

save alpha alpha